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Abstract 

There  exist  two  types  of  stationary  instability  of  the  flow  over  a  rotating  disc 
corresponding  to  the  upper-branch,  inviscid  mode  and  the  lower-branch  mode,  which  has  a 
triple  deck  structure,  of  the  neutral  stability  curve.  A  theoretical  investigation  of  the  linear 
problem  and  an  account  of  the  weakly  nonlinear  properties  of  the  lower  branch-modes  have 
been  undertaken  by  Hall  (1986)  and  MacKerrell  (1987)  respectively.  Motivated  by  recent 
reports  of  experimental  sightings  of  the  lower-branch  mode  and  an  examination  of  the 
role  of  suction  on  the  linear  stability  properties  of  the  flow  here  we  investigate  the  effects 
of  suction  on  the  nonlinear  disturbance  described  by  MacKerrell  (1987).  The  additional 
analysis  required  in  order  to  incorporate  suction  is  relatively  straightforward  and  enables 
us  to  derive  an  amplitude  equation  which  describes  the  evolution  of  the  mode.  For  each 
value  of  the  suction  a  threshold  value  of  the  disturbance  amplitude  is  obtained;  modes  of 
size  greater  than  this  threshold  grow  without  limit  as  they  develop  away  from  the  point  of 
neutral  stability. 


^  This  research  was  supported  in  part  by  NASA  Contract  No.  NASl-18605  while  the  authors  were  in  residence 
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1.  Introduction 


The  stability  of  rotating  disc  flows  has  been  the  subject  of  many  studies,  one  reason  for 
this  being  that  the  basic  boundary  layer  flow  given  by  Von  Karman’s  (1921)  exact  solution 
is  fully  three-dimensional.  Interest  also  stems  from  the  observed  similarities  in  the  type  of 
instability  occurring  in  rotating  disc  flow  work  (Gregory,  Stuart  &  Walker  1955)  and  near 
the  leading  edge  of  a  swept  back  wing,  Poll  (1978).  This  instability,  known  as  cross-flow 
instability,  is  due  to  the  inflectional  character  of  the  basic  velocity  profile  and  was  first 
studied  extensively  by  Gregory  et  al.  (1955).  They  showed  with  china  clay  techniques  that 
this  instability  took  the  form  of  a  regulaxly  spaced  pattern  of  equi-angular  spiral  vortices 
which  was  stationary  relative  to  the  disc.  Stuart  (in  Gregory  et  al.),  using  inviscid  theory, 
suggested  that  the  instability  could  be  associated  with  a  particular  inflectional  profile  in 
which  the  inflection  point  coincided  with  a  point  of  zero  velocity  somewhere  in  the  flow.  His 
calculation  gave  the  predicted  number  of  vortices  to  be  approximately  four  times  greater 
than  the  observed  value  of  about  30  but  the  angle  of  13°  between  the  axes  of  the  vortices 
and  the  radius  vector  was  in  excellent  agreement  with  their  experiments. 

Several  attempts  have  been  made  to  explain  the  observed  spiral  patterns  by  means 
of  a  linear  stability  theory  (see,  for  example,  Cebeci  k,  Stewartson  (1980),  Malik  et  al. 
(1981),  Mack  (1985),  Malik  (1986)).  Malik  et  al.  (1981)  showed  that  Coriolis  force  and 
streamline  curvature  cannot  be  neglected-  a  point  first  made  by  Faller  k  Kaylor  (1966)  and 
Lilly  (1966)  in  the  context  of  the  Ekman  boundary  layer.  The  theoretical  results  of  Mack 
(1985)  strongly  suggest  that  the  spiral  streaks  observed  in  flow  visualisation  experiments 
are  the  constant  phase  lines  of  the  merged  wave  patterns  produced  by  several  random 
sources  on  the  disc.  This  was  first  found  by  the  experimental  investigations  of  Wilkinson 
k  Malik  (1985)  who  revealed  that  the  wave  patterns  from  each  point  source  spread  out 
circumferentially  downstream  of  the  source  and  that  the  wave  patterns  from  a  number  of 
sources  eventu8dly  merge  and  cover  the  entire  circumference  of  the  disc. 

Experiments  by  Fedorov  et  al.  (1976)  showed  that  for  various  rotation  rates  of  the  disc 


an  instability  occurred  which  also  appeeired  as  a  pattern  of  spiral  vortices  which  numbered 
in  the  region  14  to  16  and  which  had  their  axes  inclined  at  angles  of  about  20°  to  the 
radius  vector.  These  discrepancies  between  these  observations  and  those  of  Gregory  et  al. 
(1955)  imply  that  the  vortices  seen  by  Fedorov  et  al.  were  not  those  of  Gregory  et  al. 
Malik  (1986)  calculated  the  neutral  stability  curve  for  stationary  disturbances  and 


demonstrated  that  two  types  of  mode  could  exist.  The  first,  the  upper  branch  of  the 
neutral  curve,  corresponds  to  the  inviscid  mode  described  by  Stuart  in  Gregory  et  al.  The 
second,  the  lower  branch  mode,  is  an  essentially  viscous  disturbance  which  corresponds  to 
zero  wall  shear  stress  of  the  crossflow  velocity  profile. 


C.-'das 
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A  linear  asymptotic  investigation  of  the  inviscid  mode  at  high  Reynolds  number  was 
conducted  by  Hall  (1986);  work  that  has  been  extended  to  include  nonlinear  effects  by 
Gajjar  (personal  communication).  Further,  Hall  (1986)  demonstrated  that  the  lower 
branch  viscous  mode  is  governed  by  a  triple  deck  structure  (which  is  derived  from  the 
classical  structure  as  reviewed  by,  for  example.  Smith  1982)  and  it  is  this  mode  which 
is  thought  to  have  been  observed  by  Fedorov  et  al.  (1976).  MacKerrell  (1987),  hereafter 
referred  to  as  M,  extended  the  study  of  Hall  (1986)  in  order  to  give  an  account  of  the  weakly 
nonlinear  properties  of  the  lower  branch  mode.  Following  the  framework  set  up  in  Hall 
(1986)  M  derived  solutions  which  enabled  the  finite  amplitude  growth  of  a  disturbance  close 
to  neutral  to  be  described.  Importantly,  she  found  that  nonlinearity  has  a  destabilising 
effect  on  the  vortices  and  discovered  the  existence  of  a  threshold  amplitude  (in  the  sense 
of  Stuart  1960,  1971  and  Watson  1960)  such  that  disturbances  smaller  than  the  threshold 
value  decay  as  they  move  away  from  the  neutral  position  whereas  those  above  threshold 
grow  explosively.  This  led  M  to  conclude  that  for  small  disturbances  the  inviscid  mode 
of  Gregory  (  '  al.  (1955)  is  probably  dominant  whereas  for  larger  amplitude  perturbations 
the  short  wavelength  mode  observed  by  Fedorov  et  al.  (1976)  is  the  more  important. 

Until  very  recently  there  have  been,  to  the  best  of  the  authors’  knowledge,  few,  if 
any,  other  observations  of  the  short  wavelength  disturbance.  At  the  1991  ICASE/LaRC 
workshop  on  transition  and  turbulence  Corke  reported  that  in  his  experiments  on  rotating 
disc  flows  a  subcritical  stationary  instability  had  been  observed-  the  mode  which  is  almost 
certainly  that  seen  by  Fedorov  et  al.  (1976)  and  analysed  in  M.  At  the  same  meeting 
we  learned  of  an  investigation  by  Dhanak  (1991)  who  has  extended  the  calculations  of 
Malik  (1986)  to  include  the  effects  of  suction.  The  role  of  suction  in  practical  flows  is  an 
important  one  for  it  is  believed  that  imposition  of  suitable  suction  on  a  boundary  layer 
can  often  stabilise  the  flow  and  thence  delay  the  onset  of  transition.  Dhanak  (1991)  found 
that  if  suction  is  applied  to  the  rotating  disc  flow  then  the  overall  effect,  at  least  according 
to  linear  theory,  Is  that  the  flow  is  indeed  stabilised:  conversely  blowing  destabilises  the 
flow.  This  finding  invites  the  question  as  to  the  effect  of  suction  on  the  nature  of  the  short 
wavelength  mode  described  by  Hall  (1986)  and  M  and  it  is  this  that  we  address  here. 

In  §2  we  follow  the  nonlinear  analysis  of  M  and  indicate  the  modifications  required  in 
order  to  incorporate  suction  into  the  flow.  We  then  execute  some  numerical  work  to  obtain 
a  quantitative  description  of  the  flow  properties  and  conclude  with  some  discussion. 
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2.  Formulation  of  the  problem 

We  consider  the  case  in  which  the  disc  z  =  0  rotates  about  the  z-axis  with  angular 
velocity  Q.  Relative  to  cylindrical  polar  co-ordinates  [r,6,z)  which  rotate  with  the  disc, 
the  continuity  and  Navier-Stokes  equations,  when  suitably  non-dimensionalised,  become 

V.u  =  0,  (2.1) 

^  +  (u. V)u  +  2(k  A  u)  -  rf  =  -  Vp  +  —  V^u,  (2.2) 

ot  lie 

where  u  =  {u,v,w)  sire  the  velocity  components,  p  is  the  fluid  pressure,  f  and  k  are  unit 
vectors  in  the  ] —  and  z—  co-ordinate  directions  and  the  Reynolds  number  Re  =  /u 
with  u  the  kinematic  viscosity  of  the  fluid  and  L  a  reference  lengthscale.  The  Reynolds 
number  is  taken  to  be  large  throughout. 

For  the  ensuing  analysis  it  is  convenient  to  define  e  =  Re~^  and  with  z  =  e^r}  the 
basic  steady  flow  is  given  by  the  solution 

u  =  UB  =  {ru(r)),  rv{r)),  £^w(t]))  ,  p  =  pg  ^  p{r)), 
where  u,  v,  w  and  p  satisfy 

—  (1  +  v)'^  +  u'w  =  u'\  2u(l  +  u)  +  v'w  =  v” ,  (2.3a,  b) 

w'  +  2u  =  0,  p'  +  ww'  —  w"  =  0.  (2.3c,  d) 

The  appropriate  far  field  conditions  take  the  forms 

u  — ►  0,  V  — *  —1  as  Tj  — >  oo,  (2.4a) 

whereas  on  the  disc  surface  we  impose 


u  =  u  =  0, 


(2.46) 


together  with 

w  =  —S,  (2.4c) 

where  S’  is  a  measure  of  the  suction  (S  >  0)  or  blowing  (S  <  0)  applied  to  the  flow. 

The  linear  stability  of  the  basic  flow  (2.3),  (2.4)  in  the  absence  of  suction  or  blowing 
was  addressed  by  Hall  (1986).  He  demonstrated  that  the  fiov  is  susceptible  to  a  stationary 
disturbance  which  is  governed  by  an  appropriate  triple  deck  structure  whose  lower,  main 
and  upper  decks  are  of  thicknesses  0(£®),  0(£*)  an  1  0{£*)  respectively.  Hall  showed  that 
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the  disturbance  wavenumbers  in  the  r  and  9  directions  axe  0{e  and  solutions  were 
sought  proportional  to 

E  =  exp  a{r,e)dr+ 9 /3{e)^  ,  (2.5) 

where  the  wavenumbers  a  and  j3  expand  as 

Q  =  &0  +  +  £^Q2  +  .  .  .  ,  P  =  Pq  Pi  +  P2  +  .  .  .  .  (2-6) 


It  was  found  by  Hall  (1986)  that,  as  shown  by  Gregory  et  al.  (1955),  the  ‘effective’ 
velocity  profile  for  a  three-dimensional  disturbance  with  wavenumbers  a  and  P  in  the  r 
and  9  directions  is  rau  +  Pv.  However,  if  the  effective  wall  shear  rau'  +  Pv'  does  not  vanish 
the  modes  are  necessarily  time  dependent  for  a  and  P  real.  Therefore,  in  order  to  seek 
stationary  modes  which  are  neutrally  stable  at  the  location  r  the  values  01  j  and  Pj  must 
be  chosen  so  as  to  reduce  the  effective  wall  shear  to  zero  at  leading  order. 

MacKerrell  (1987)  extended  the  work  of  Hall  (1986)  to  a  weakly  nonlinear  setting.  She 
demonstrated  that  in  order  to  obtain  a  classical  evolution  equation  for  the  disturbance  it 
is  necessary  to  examine  disturbances  of  amplitude  0{£  5 )  within  the  lower  deck.  Further, 
the  mode  then  evolves  in  an  0(e)  neighbourhood  of  the  neutral  position  r  =  f  and,  as  in 
usual  nonlinear  triple  deck  work,  it  is  the  lower  tier  of  the  triple  deck  where  nonlinearity 
manifests  itself.  Our  objective  here  is  to  generalise  the  work  of  M  to  investigate  the  role 
of  the  suction  parameter  S  and  in  the  following  we  deliberately  give  the  briefest  of  details. 
The  analysis  differs  only  slightly  from  that  of  M  and  so  we  highlight  the  main  differences 
between  the  present  work  and  that  contained  in  M.  The  omitted  details  may  be  found  in 
M  and  the  reader  is  referred  to  that  paper  where  necessary. 

Following  the  argument  of  M  we  define  a  new  variable  ri  in  the  neighbourhood  of  the 
neutral  point  f  by 

r  =  f  +  eri,  (2.7a) 

and  seek  disturbances  whose  radial  dependence  may  be  written  in  the  form 

d  iao  iai  1  d 
dr  e  dri 

Within  the  lower  deck,  where  we  define  the  0(1)  co-ordinate  ^  by  2  =  we  suppose 
that  the  basic  flow  velocity  components  (2.3)  expand  according  to 


u  =£Uo^  +  +  ■■■, 

V  =£Vo^  +  e^vi^^  +  +  . . . ,  (2.8) 

u;  =  -  5  +  +  e^u)2^^  -f - 
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If  we  perturb  the  basic  flow  by  writing 


{u,v,w,p)  =  (u  +  U,v  +  V,e^w  P), 

then  in  the  lower  deck  the  disturbance  quantities  develop  according  to 


U  =£j 


(^0  +  4-  •  .  •)  "I - h  +  •  •  •  E 

s/3S  £  J 

+  harmonic  terms  +  complex  conjugate, 


together  with  a  similar  expansion  for  V.  In  addition 


W  =e^ 


7  \  noe^Ci 


[(oto^ui  +  /doUl  +  . . .]  +  e^Wi  +  W2  +  . . .  E 
+  harmonic  terms  +  complex  conjugate, 


(2.10a) 


(2.106) 


P=£?  [{e^Pi+e*P2  +  ...)]E 

+  harmonic  terms  +  complex  conjugate. 


(2.10c) 


Here  Ci(ri)  is  a  scaled  amplitude  of  the  disturbance  and  we  have  made  the  definition 


-2  _  2  ,  ^0 

7o  =  Oq  + 


(2.11) 


These  forms  of  the  disturbance,  in  which  all  the  unknowns  are  functions  of  ^  and  ri ,  were 
first  derived  by  M  who  also  elucidated  the  forms  of  the  harmonic  terms.  However,  we  find 
that  the  crucial  harmonic  and  mean  flow  terms  which  interact  to  drive  the  fundamental 
and  hence  lead  to  the  desired  amplitude  equation  are  unchanged  by  the  imposition  of 
suction  or  blowing.  Hence,  for  brevity  we  do  not  repeat  these  results  from  M  and  now 
merely  state  how  her  overall  findings  need  to  be  modified. 

Firstly,  following  the  procedure  adopted  in  M,  we  find  that  W2  satisfies 

- s^-g^+2sW2-A  -2ri~f^Ci-^ +  i-fQVol3iA  * 


—  i'?  rt  — i.  vjo 

+  iki{faiUo  +  0ivo)^  *£  ~ip2VoA 


-llC2Ct- 


-  2i7oriCi 


0lr 

■  17(0,  V2») 

,  alJ  C7(0,0) 


__  A  -4  r,o\\7oDl  ,  Dtrc 

-  ruoA  — +2SA*—^  +  RHS, 

00  dri  01 


(2.12a) 
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where  A  =  i{aofui  +  PqVi),  Di{ri)  and  C2(ri)  are  functions  of  Ci  given  in  M,  ki  is  a  real 
constant  whose  precise  value  will  turn  out  to  be  unimportant,  s  =  and  an  asterisk  on 
a  quantity  denotes  the  complex  conjugate  of  that  quantity.  Further,  U{a,x)  denotes  the 
usual  parabolic  cylinder  function  (see,  for  example,  Abramowitz  &  Stegun  (1964))  and 


- 


RHS  =  —  2ir-^UoA  *—ikiA  ^  s^(raoU2  +  0oV2)  +  2i'yoCi—;z — -suofo 

Po  Po 

-  i(faiUo  +  PiUo)A“t  +  3iA~^s^(raoU2  +  PoU2)j 


A-^  _  _  U{0,y/2s) 


X  ^,^CiFi(s)  +  -^Ci  1  +  4 


^^0^  UoF2(s) 


Po  ^  V  ^IJ  ^(0,0) 


+  i(rooW2  +  Poi'2)A  +  i{raiUo  + 


-5A-' 


(i  +  ?l 

V  «0 


tio4'(5) 


IJ  U{0,0)\' 


where  Fi(s)  and  F2(s)  satisfy 


Fi"  -  s^F[  +  2sFi  =  1,  Fi(0)  =  Fi(oo)  =  0, 


F(0,0) 


(2.126) 


4"  -  +  2sF2  =  F(0,  sy/2),  F2(0)  =  ^2(00)  =  0. 


The  solution  for  Fo  when  S  =  0,  say  Foo?  is  given  in  Hall  (1986).  It  is  then  straightforward 
to  obtain 

r-r  7-r  ,  o  a -4  ^  ^  7o ^0  •sF(0,  S>/2) 

C/,=C;.o+SA  c;(o,o)  • 

In  order  for  the  disturbance  to  vanish  on  the  disc  surface  we  find  from  the  continuity 
equation  that  dW2lds  =  0  on  s  =  0.  By  solving  (2.12)  subject  to  this  condition  we  obtain 
the  desired  amplitude  equation  in  the  form 


—  (a  +  ic)riC\  +  (6  +  id)C\  |Cj  |^  +  (e  +  if)C\ , 


(2.13) 


where  a,  6,  c,  d,  e  and  /  are  real  valued  constants,  the  more  important  of  which  are  detailed 
presently.  Multiplying  (2.13)  through  by  C*,  adding  the  complex  conjugate  and  replacing 
ri  by  r'j  =  ri  +  (e/a)  yields  the  amplitude  evolution  equation 


-(|C,H  =  2ar;|C,P+26|C,|‘, 


(2.14) 
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and  the  form  of  the  solution  of  (2.14)  depends  on  the  signs  of  a  and  b.  ^From  (2.3)  and 
(2.4)  we  find  that 

1  Suo  ,  _  Svo 

=  and  vi  = - - 

and  recalling  the  imposition  that  the  leading  order  wall  shear  (=  aofuo  +  /douo)  must  be 
necessarily  zero  in  order  to  obtain  stationary  vortices  we  have  A  =  —ifao/2.  Following 
the  manipulations  of  M  we  finally  obtain 


f‘^uo(aof)'! 


(2.15o) 


and 


where 


,  _  3(o'or)'7o  /  h 


(2.156) 


I 


1 


and 


47o/?o-^i  -  (^)  (roiuo  + /?iUo)|  +  To^o-^i  -  (^)  (roiuo  +  ^iuo) 

h=F[iO)=  f 

Jo 

r.  -  2^2(0)  _  f 

Jo 


(2.15c) 

(2.15d) 

(2.15e) 

U{0,0) 

Consequently,  we  now  have  the  outcome  that  the  formulae  (2.15)  for  the  coefficients  of 
the  amplitude  equation  appear  to  be  independent  of  the  suction  parameter  5  and  hence 
are  identical  to  those  of  M.  However,  the  constants  appearing  in  the  expansions  (2.8)  for 
the  basic  flow  quantities  close  to  ^  =  0  are  fujictions  of  S  and  thence  so  are  a  and  6.  In 
the  following  section  we  conduct  a  few  calculations  for  various  values  of  S  and  draw  some 
conclusions. 


3.  Results  and  discussion 


In  order  to  calculate  the  coefficients  a  and  6  in  (2.14)  we  merely  need  to  solve  the 
basic  equations  (2.3)  subject  to  boundary  conditions  (2.4)  and  note  that  the  function 
(foiuo  +  di^o)  appearing  in  definition  (2.15c)  is  given  by 


To-^i  + 


tjoup 

^0 


,-.2 


1  + 


Ur 


I  _  _  _  7o 

I2  =  iA5  (otiruo  +  diUo) 

yo 
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a  relation  which  was  deduced  in  M.  This  relation  also  serves  to  determine  the  value  of  70  and 
this,  in  addition  to  the  requirement  that  the  leading  order  wall  shear  afuo  +  /Jo^o  vanishes, 
furnishes  the  leading  order  radial  and  azimuthal  wavenumbers  Oq,  Pq.  The  equations  for 
the  basic  flow  were  solved  using  a  suitable  NAG  routine  and  then  we  found  it  convenient 
to  write 

a  =  Af~^,  b  =  Bf~^,  (3.1) 


where  the  functions  A(5)  and  B{S)  are  given  in  Figure  (1).  We  note  that  for  all  values  of 
suction  or  blowing,  —  00  <  5  <  00,  we  have  A  <  0  and  B  >  0  and  the  former  condition 
implies  that  the  amplitude  of  the  solution  increases  or  decreases  depending  on  whether  r 
is  less  than  or  greater  than  the  neutral  value.  It  is  also  straightforward  to  demonstrate 
that  if  we  write 

|Cil^  =  a/26^  y  and  r[  =  xj\/—a 


then  (2.14)  becomes 


whose  solution  is 


dx 


-2xy  + 


y{x)  = 


2 


where  j/o  =  1/(0)  and  erf(a;)  =  dt.  Since  y{x)  must  be  greater  than  or  equal  to 

zero  and  0  <  erf(x)  <  1,  the  solution  is  bounded  and  valid  for  all  x  if  yo  <  2/y^.  However, 
for  yo  >  2/v/ff,  y(x)  becomes  infinite  at  some  point  xq,  0  <  xq  <  00,  and  the  solution  is 
only  valid  for  x  <  xq.  Thus  the  value  yo  =  2/^/^  represents  a  threshold  between  solutions 
which  grow  indefinitely  and  those  that  decay  to  zero.  The  corresponding  disturbance 
threshold  value 


(3.2) 


is  shown  in  Figure  (2).  We  note  that  for  large  suction  (5  »  1)  the  threshold  amplitude 
becomes  small  whereas  for  large  blowing  this  amplitude  grows.  We  now  verify  these 
qualitative  results  by  considering  the  appropriate  limits. 


3.1.  The  large  suction  limit,  5  >  1. 


For  5  1,  Stuart  (1934)  has  demonstrated  that  the  basic  flow  profile  is  confined 

within  a  zone  of  thickness  0(1/5)  and  if  we  write  ^  =  517  where  ((  =  0(1)  then,  to  leading 
order, 


u  = 


252 


,-c 


)  ,  V  =  1  —  e 


and  w  =  —S. 
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Consequently  the  constants  defined  in  the  basic  flow  expansions  (2.8)  assume  the  forms 

1  -  ^ 

uo  =  ^0  =  -S 

and  these  lead  to  the  asymptotic  expressions 

A  ~  «  -1.20785^,  (3.3a) 

X 

24/5  -  27  -  27 

B - =311 - g-5V  w  8.19525^,  (3.36) 

(21  -  10\/2)J/ 

where  Ii  and  are  given  by  (2.136,  c)  and  A  and  B  are  the  coefficients  in  the  amplitude 
equations  (2.14)  and  which  were  defined  in  (3.1).  In  turn,  the  threshold  amplitude 

|Ci|e  w  0.3063S~5f5.  (3.3c) 


3.2.  The  large  blowing  limit,  S  <C  —1. 

In  his  consideration  of  boundary  layer  flows  Pretsch  (1944)  demonstrated  that  the 
influence  of  viscosity  becomes  almost  negligible  for  sufficiently  large  blowing  and  this  work, 
together  with  the  results  emanating  from  studies  of  fluid  injection  in  supersonic  boundary 
layers  by  Smith  &  Stewartson  (.1973a,  6),  enable  us  to  derive  the  flow  structure  for  the 
present  problem  when  5  <C  —  1.  We  find  it  convenient  to  define  the  positive  quantitv 
S  —  —S.  It  is  now  the  case  that  the  extent  of  the  boundary  layer  increases  and  formally 
becomes  0{S).  Therefore,  if  we  define  the  0(1)  co-ordinate  Zi  by  2i  —  rj/S  we  find  that 
in  the  region  0  <  2i  <  7r/2  the  basic  flow  solution  develops  according  to 


U  —  Fo(^i)-|-^Fi(2i)  +  ...,  U  —  Go(^i)+^Gi(2i)  +  ..., 

5  S 


(3.4) 


with 


Hq  =  Gq  =  COS^ 


and  Fq  =  cos  zi  sin  zi . 


(3.5) 


The  profile  decays  algebraically  as  zi  — ►  it  12  and  a  thin  viscous  layer  centred  on  this 
point  ensures  that  the  solutions  for  u  and  v  are  brought  to  zero  at  infinity.  Details  of  this 
viscous  structure  are  not  needed  for  present  purposes  as  we  are  solely  concerned  with  the 
values  of  u'(0)  and  v'(0).  Consideration  of  the  problem  for  Gi(2i)  enables  us  to  conclude 
that  the  flow  constants  uq  and  uo  in  (2.8a,  6)  are  given,  at  leading  orders,  by 


and 
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In  turn  this  leads  to  the  asymptotic  results 

^4  ^  «  -0.3020  (-5)^  ,  (3.6a) 

i 

6/4  _  _ii 

B - ^ - tS  "  «  2.0488  (-S)  ^  ,  (3.6i) 

(21  -  1QV2]I2* 

with  a  threshold  amplitude  given  by 

|Ci|e  0.3063(-5)^  ft.  (3.6c) 

4.  Conclusions 

On  Figures  (1)  &  (2)  we  have  indicated  the  asymptotic  behaviours  (3.3)  &  (3.6)  both 
for  5  <  —1  and  for  5  >  1.  We  observe  that  these  one  term  asymptotic  predictions  are 
very  accurate  for  a  surprisingly  large  range  of  values  of  5;  indeed  they  are  graphically 
indistinguishable  from  the  respective  curves  for  |5j  greater  than  about  five. 

The  most  important  deduction  to  be  drawn  from  Figure  (16)  is  that  the  nonlinear 
coefficient  in  evolution  equation  (2.14)  is  positive  for  all  values  of  blowing  and  suction. 
Thus  in  all  cases  nonlinearity  has  a  destabilising  influence  and  we  have  found  the  existence 
of  a  threshold  amplitude  |Ci|e  as  discussed  by  Stuart  (1960,  i971)  and  Watson  (1960)  such 
that  all  disturbances  of  amplitude  less  than  threshold  decay  whereas  those  greater  than 
threshold  become  infinite  as  the  distance  from  the  neutral  stability  position  increases.  This 
in  turn  leads  to  turbulence. 

Our  results  described  above  seem  to  suggest,  perhaps  unexpectedly,  that  the  effect  of 
including  suction  in  the  rotating  disc  problem  tends  to  lower  the  threshold  amplitude  for 
the  disturbance,  .«ee  Figure  (2)  and  result  (3.3c).  Moreover,  Figure  (2)  indicates  that  only 
very  moderate  suction  or  blowing  is  needed  in  order  to  have  a  significant  effect  on  the 
critical  disturbance  amplitude.  In  the  work  of  Dhanak  (1991)  it  was  found  that  suction 
stabilises  the  flow  in  as  much  thnt  inclusion  of  suction  reduces  the  region  of  wavenumber/ 
Reynolds  number  parameter  space  in  which  the  basic  profile  is  unstable  to  stationary, 
infinitesimal  modes.  However,  since  the  threshold  amplitude  decreases  with  increasing 
suction  we  speculate  that  an  experimental  configuration  would  need  to  be  less  strongly 
forced  in  order  for  the  subcritical  instability  to  occur  than  would  be  the  case  for  the  zero 
suction  problem:-  conversely,  imposition  of  blowing  destabilises  the  flow  according  to  linear 
analysis  but  the  much  larger  threshold  amplitude  suggests  that  exciting  the  subcritical 
disturbance  experimentally  might  be  extremely  difficult. 
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We  now  say  a  few  words  about  the  extension  of  our  present  problem  into  the 
compressible  regime.  Seddougui  (1990)  considered  the  effects  of  compressibility  on  the 
subcritical  mode  of  M  and  found  that  for  both  an  adiabatic  wall  and  an  isothermal  wall 
the  stationary  mode  is  only  possible  over  a  finite  range  of  Mach  numbers.  The  analysis 
conducted  by  Seddougui  (1990)  is  naturally  more  involved  than  that  of  M  but  we  have  seen 
here  that,  at  least  for  the  incompressible  flow  problem,  the  inclusion  of  suction  requires 
only  relatively  minor  changes  to  the  nonlinear  analysis.  With  suction  the  variations  in 
the  coefficients  of  amplitude  equation  (2.14)  axe  predominantly  due  to  the  changes  in  the 
wall  shears;  these  are  in  turn  due  to  the  changing  boundary  condition  (2.4c)  for  the  basic 
flow.  Preliminary  consideration  of  the  compressible  version  of  the  present  problem  also 
suggests  that  the  modifications  required  in  order  to  account  for  suction  in  Seddougui ’s 
(1990)  work  will  also  be  relatively  straightforward.  Further  studies  on  this  aspect  are  in 
progress  together  with  an  investigation  into  the  properties  of  the  time  dependent  version 
of  the  mode  considered  here. 

Finally  there  is  the  question  of  the  extension  of  the  present  work  into  the  fully  nonlinear 
regime.  Examples  of  fully  nonlinear  calculations  for  suction/blowing  problems  are  typified 
by  those  of  Smith  &  Stewartson  (1973a,  6)  for  plate  and  slot  injections  into  supersonic 
boundary  layers.  Extension  of  our  work  along  the  lines  suggested  in  these  papers  would 
be  valuable  indeed. 
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^ogioi-A) 


Figure  (1).  The  dependences  of  amplitude  equation  coefficients  A  EUid  B  defined  by  (2.14), 
(2.15)  and  (3.1)  on  the  suction  parameter  S.  a)  logio(— >1):  dotted  lines  indicate  the 
asymptotic  values  (3.3a)  and  (3.6a),  logjo^:  dotted  lines  indicate  the  asymptotic  values 
(3.36)  and  (3.66). 
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